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NARROW SIEVES FOR 
PARAMETERIZED PATHS AND PACKINGS 

ANDREAS BJORKLUND 1 , THORE HUSFELDT 2 , PETTERI KASKI 3 , 
AND MIKKO KOIVISTO 4 

Abstract. We present randomized algorithms for some well-studied, 
hard combinatorial problems: the fc-path problem, the p-packing of q- 
sets problem, and the g-dimensional p-matching problem. Our algo- 
rithms solve these problems with high probability in time exponential 
only in the parameter (fc, p, q) and using polynomial space; the con- 
^^ ■ stant bases of the exponentials are significantly smaller than in previous 

C^ works. For example, for the fc-path problem the improvement is from 2 

to 1.66. We also show how to detect if a d-regular graph admits an edge 
coloring with d colors in time within a polynomial factor of 0(2^ d ~ 1 ' n ' 2 ). 
Our techniques build upon and generalize some recently published 
O ideas by I. Koutis (ICALP 2009), R. Williams (IPL 2009), and A. Bjork- 

lund (STACS 2010, FOCS 2010). 

q 

q ■ 1. Introduction 

Combinatorial problems such as finding a long simple path in a graph or 
disjointly packing many members of a set of subsets are well-studied and 
hard. In fact, under standard complexity-theoretic assumptions, algorithms 
\^Q \ for these problems must either be inexact, or require running times of super- 

polynomial or maybe even exponential time. 

It has been observed that the complexity of the problems studied in the 
present paper depends exponentially on the output size k instead of the input 
size n, i.e., they admit running times of the form exp(poly(A:))-poly(n) rather 
than, say, exp(0(n)). A number of papers have improved the exponential 
dependencies dramatically over the past decade, arriving at exponential fac- 
tors of size 2 . We further improve this dependency, reducing the exponent 
rS | base below the constant 2, sometimes significantly. 

To express our results in the terminology of parameterized computational 
complexity theory, we improve the running time of several canonical fixed 
parameter tractable problems. In particular, we claim the ephemeral lead in 
the highly competitive "FPT races" for path finding, uniform set packing, 
and multidimensional matching. 

Our techniques also allow us to report progress on an unparameterized 
problem: we show a nontrivial upper bound on the complexity of edge col- 
oring for regular graphs. 

We adopt the notational convention from parameterized and exponential 
time algorithms, letting 0*(f(k)) denote f{k)n°^ l \ where typically n is 
some aspect of the input size such as number of vertices, and A; is a parameter 
such as path length. Our parameters are all polynomially bounded by n, 
so O* also hides factors that are polynomial in the parameter size. We 
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present our results in terms of decision problems, they can be turned into 
optimization or search problems by self-reductions in the obvious way. 

All our algorithms are randomized. The error is one-sided in the sense that 
they never report a false positive. The error probability is constant and can 
be made exponentially small by a polynomial number of repetitions, which 
would again be hidden in the O* notation. 

1.1. Finding a path. Given an undirected graph G on n vertices, the k- 
path problem asks whether G contains a simple path on k vertices. 

Theorem 1. The undirected k-path problem can be solved in time 0*(1.66 fc ) 
by a randomized algorithm with constant, one-sided error. 

The proof is in $2j For k = n, the result matches the running time of the 
recent algorithm for Hamiltonian path of Bjorklund [3j. 

Previous work. Naively, the fc-path problem can be solved in time 0*(n k ), 
but Monien [29] and Bodlaender [5] showed that the problem can be solved 
in time 0*(f(k)), leading Papadimitriou and Yannakakis [30J to conjecture 
that the problem was polynomial-time solvable for k = O(logn). This was 
confirmed in a strong sense by Alon, Yuster, and Zwick, with a beautiful 
0*(c ) algorithm. A number of paper have since reduced the base c of the 
exponent using different techniques, see Table 1. (In the following tables, 
we mark randomized algorithms with 'r'.) 



Table 1. 


fc-path in time 0*(f(k)) 




k\ 


Monien [29] 


12.6 fe Chen et al. [6] 


m k 


Bodlaender [5] 


4 fe r Chen et al. [B] 


5.44 fc r 


Alon et al. Q] 


2.83 fc r Koutis (2008) gj 


c k 


c > 8000, Alon et al. pQ 


2 k r Williams [37J 


16 fc 


Kneis et al. [25J 


1.66 fc r this paper 



Our result is yet another improvement of the exponent base, notable per- 
haps mostly because it breaks the psychological barrier of c = 2. We fully 
expect this development to continue. On the other hand, computational 
complexity informs us that an even more ambitious goal may be quixotic: 
An algorithm for /c-path with running time exp(o(/c)) would solve the Hamil- 
tonian path problem in time exp(o(n)), which is known to contradict the 
exponential time hypothesis [18] . 



1.2. Packing disjoint triples. Let 9" be a family of subsets of an n-element 
ground set. A subset A C 3" is a p-packing if \A\ = p and the sets in A are 
pairwise disjoint. Given input set 3~ of size-3 subsets, the p-packing of 3- 
sets problem asks whether 3" contains a p-packing. This problem includes 
a number of well-studied problems in which the ground set consists of the 
vertices of an input graph G = (V,E). In the vertex-disjoint triangle p- 
packing problem, the set 3" consists of the subsets of V that form a triangle 
K3. In the edge-disjoint triangle p-packing problem, the set 3" consists of 
the subsets of E that form the edges of a triangle. In the vertex-disjoint P3 
p-packing problem, the set J consists of the vertex subsets {u, v, w} C V for 
which uv, vw G E. 



Theorem 2. The p-packing of3-sets problem can be solved by a randomized 
algorithm in time 0*(1.493 3p ) with constant one-sided error. 

The proof is in £0] For p = n, the result matches the running time of the 
recent algorithm for exact cover by 3-sets of Bjorklund [2J. 

Previous work. The naive algorithm for p-packing considers all (' ') ways of 
selecting p sets from 3~. Results of the form 0*(f(p)) go back to Downey 
and Fellows [7], and the dependency on p has been improved dramatically in 
a series of papers, see Table 2. Remarkably, the best previous running time 
is given by Koutis's algorithm for the more general problem of p-packing 
sets of size q, specialized to the case q = 3. (We return to the performance 
of our own algorithm in the case q > 3 in §1.31 ) 



Table 2: p 


-packin; 


gs of 3-sets in time 0*(f(p)) 


2°W(3p)! 






Downey and Fellows [7] 


(5.7p) p 




r 


Jia, Zhang, and Chen [19] 


2°(p) 






Koutis (2005) [2DJ 


(12.7c) 3p 






c > 10.4, Fellows, Heggernes, et al. [TTJ 


10.88 3p 




r 


Koutis (2005) [2DJ 


4.68 3p 






Kneis et al. [25] 


^3p+o(p) 






Chen et al. [6] 


2.52 3 p 




r 


Chen et al. [pj 


22plogp+1.869p 




vertex-disjoint triangles K3, Fellows, Heggernes, et al. [TTJ "1" 


22.628 plogp+p 




edge-disjoint triangles K3, Mathieson et al. [28 


4.61 3p 






Liu et at. [27] 


3.523 3p 






Wang and Feng [35] 


3.404 3 p 






vertex-disjoint paths P3, Prieto and Sloper [2JJ 


2.604 3 p 






vertex-disjoint paths P3, Wang et al. [3DJ 


2.482 3 ^ 






vertex-disjoint paths P3, Fernau and Raible [T5] 


2 3p 




r 


Koutis [2TJ 


1.493 3p 




r 


this paper 



t The precise time bound can be seen to be 0(2 2 P 1 °sP+ 1 - 869 P n 3 ). 

It is known that packing vertex-disjoint copies of H into G is NP-complete 
as soon as H is connected and has more than 2 vertices [15] . Fellows et al. 
[11] raised the question of how hard the parameterized problem is, and we 
observe here that this can be answered under the exponential time hypoth- 
esis [18] . which is equivalent to the parameterized complexity hypothesis 
FPT^Mfl]. 

Proposition 3. There is no algorithm for vertex-disjoint triangle p-packing 
in time exp(o(»)) unless the satisfiability of 3- CNF formulas in n variables 
can be decided in time exp(o(n)). 

The proof of the proposition is routine and we omit a detailed presen- 
tation. Briefly, a 3-CNF instance to satisfiability is transformed to a 3- 
dimensional matching instance; the standard reduction results in a graph 
of size 0{n + m), and the sparsification lemma of Impagliazzo, Paturi, and 
Zane [18] is used to remove the dependency on m. A subexponential time 
(in p < n) algorithm for this problem would solve the original instance in 
time exp(o(»)) = exp(o(n)). 



1.3. Uniform set packing. As before, let 3" be a set of subsets of an n- 
element ground set. A subset A C 3" is a p-packing if \A\ = p and the 
sets in A are pairwise disjoint. Given as input a family 3" of size-g subsets, 
the p-packing of q-sets problem asks us to determine whether 3" contains a 
p-packing. 

This is a generalization of the triple p-packing problem described in N1.21 
and the algorithm we advertised in that section is merely a specialization of 
a more general result. 

Theorem 4. The p-packing of q-sets problem can be solved by a randomized 
algorithm in time 0*(f(p,q)) with constant, one-sided error, where 

' 0.108157 • 29(1 - 1.64074/g) 1 - 64076 -^ - 679625 ^ p 



JyPiQ) 1 / -|\0.679623 

Potentially, 3" can have size ( n ), so reading in the input alone can take 
super-polynomial time in n. Thus we adopt the convention that |3~| is poly- 
nomial in n. Thus, the O* notation suppresses polynomial factors in n (and 
hence in p, q < n) and also in |3~|; a more careful (and even less readable) 
bound on the running time is given in £ )2.131 

Still, the above expression is difficult to parse. The previous best bound 
is Koutis's [21] much cleaner 0*(2 qp ), and our bound is not 0*((2 — e) qp ) 
for any e. Instead, our algorithm behaves well on small q; for comparison, 
we can express bounds on / of the form 0*{c qp ) for small q > 3. 

q 3 4 5 6 7 8 

f{p,q) 1.4953 3p 1.6413 4p 1.7205 5p 1.7707 6p 1.8055 7p 1.8311 8p 

q 10 20 50 100 500 

f(p,q) 1.8663 10p 1.9345 20p 1.9741 50p 1.9871 100p 1.9975 500p 

The proof is in £JH For p = n, the result matches the running time of the 
recent algorithm for exact cover by g-sets of Bjorklund [2]. 

Previous work. Most of the work on p-packing of q-sets has been done for the 
special case q = 3, described in §T72J For general q, the first algorithm with 
running time of the form 0*(f(p, q)) is due to Jia et al. [T9]; the subsequent 
improvements are shown in Table 4. 



Table 4: p-packing of g-sets in time 0*(f(p,q)) 

exp(0(pq)) Fellows, Knauer, et al. [TJ] 

5.44<?p r Koutis [2D] 

2« p r Koutis ED 



For the special graph packing case where 3" consists of isomorphic copies 
of a fixed graph H of size q, an earlier result established 0(2 p,?logp+2p9log9 n' ? ) 
|llj . The only specific packing problem we are aware of that our general 
algorithm does not seem to improve is the problem of packing vertex-disjoint 
stars into a given graph. A star Ki^-i consists of a center vertex connected 
to q — 1 other vertices. Prieto and Sloper |31j exhibit a kernel of polynomial 
size s(p, q) = p(q 3 + pq 2 + pq + 1) for this problem. They do not express 



running times in terms of q, but their "brute force" algorithm can be seen 
to run in time within a polynomial factor of ( s{ - p,q >") = exp (0(p log pq)). 

Significant improvements for p-packing of q-sets, such as a general exp(q ■ 
o(p)) algorithm, are ruled out by Proposition [3] For the nonuniform p- 
packing problem, when there is no bound on the size of the packed sets, 
we are unlikely to find an algorithm with running time 0*{f{p)) for any 
function /; the main evidence is provided in terms of parameterized com- 
plexity, where the more specific problem of finding an independent set of 
size p (equivalently, packing p subsets of is the family 3" of closed vertex 
neighborhoods in a given graph) is VK[l]-hard [7]. 

1.4. Multidimensional matching. Let U±, U2, ■ ■ ■ , U„ be pairwise disjoint 
sets, each of size r. Let 3" be a set of subsets of U\ U U2 U • • • U U q such that 
each A € 3" satisfies \A D Uj\ = 1 for each j = 1, 2, . . . , q. Given 3~ as input, 
the q-dimensional p-packing problem asks whether 3" contains a p-packing. 
One often views this problem as (/-dimensional p-matching, in which case 
3" is thought of as the edge set of a (/-uniform (/-partite hypergraph over 

JJ\ X Ui X • • • X U q . 

Again, 'J itself can have size up to r q , so reading the input alone would 
take time super-polynomial in the size n = qr of the universe already. Thus, 
we adopt the convention that |3~| is polynomial in n. In particular, the O* 
notation hides factors polynomial in both n (and hence p,q < n) and \3^\. 

Theorem 5. The q-dimensional p-packing problem can be solved in time 
0*{2^ q ~~ 2 > p ) by a randomized algorithm with constant, one-sided error. 

The proof is in Ej3j For p = r, the result matches the running time of the 
recent algorithm for (/-dimensional perfect matching of Bjorklund [2]. 

Previous work. The first parameterized multi-dimensional matching algo- 
rithm appears to be Downey, Fellows, and Koblitz's [8] application of the 
color-coding technique. Some of the ensuing improvements apply only to 
the 3-dimensional case. See Table 3. 



Table 3: (/-dimensional p-packing in time 0*(f(p,q)) 



(qp)KlP) 3qP+1 


Downey, Fellows, Koblitz [5] 


exp(0(pq)) 


Fellows, Knauer, et al. [T2] 


exp(0(pq)) 


Koutis (2005) [23 


(4e)i>s 


r Koutis (2005) [23 


2.80 3p 


q = 3, Chen et al. [B] 


2.773* 


q = 3, Liu et al. [37] 


2.52 3 ^ 


r q — 3, Chen et al. [B] 


2 qp 


r Koutis (2008) [21] 


2(<?-i)p 


r Koutis and Williams [35] 


2 (<?-2)p 


r this paper 



1.5. Edge colouring. Finally, we turn to an unparameterized problem. 

Let G be an undirected loopless n-vertex (i-regular graph without parallel 
edges. The edge-coloring problem asks whether the edges of G can be colored 



so that no two edges that share an endvertex have the same color. It is well 
known that the number of colors required is either d or d + 1. 

Theorem 6. The d-edge coloring problem for d-regular graphs can be solved 
in time O* {2^ d ~ 1 ' n ' 2 ) and polynomial space by a randomized algorithm with 
constant, one-sided error. 

The proof is in SJ3J 

Previous work. The naive algorithm for edge coloring tests all d m assign- 
ments a d colors to m edges. To the best of our knowledge, the only 
other known exact algorithm is to look for a vertex d-coloring of the line 
graph L{G) of G. The line graph L(G) has m vertices, so the algorithm of 
Bjorklund, Husfeldt, and Koivisto |4J solves the problem in time (and space) 
0*(2 m ), which is 0*(2 ' 2 ) for d-regular graphs. For 3-coloring, the current 
best bound is 0(1.344") |24j : very recently, exponential space enumeration 
algorithms for 3, 4, and 5-edge colouring were announced [14] . with running 
time O(1.201 n ), <3(1.8172 n ), and 0(3.6626 n ), respectively. 

It is known that for each d > 3 it is an NP-complete problem to decide 
whether d colors suffice |17[ l26j. However, the exponential time complexity 
of edge coloring remains wide open [23] . Curiously, we do not know how 
to apply these ideas to vertex coloring; a polynomial space algorithm with 
running time 0*(2 n ) has not yet been found. 

1.6. Methods. Our methods follow the idea introduced by Koutis [21] of 
expressing a parameterized problem in an algebraic framework by associat- 
ing multilinear monomials with the combinatorial structures we are looking 
for, ultimately arriving at a polynomial identity testing problem. Various 
ideas are used for sieving through these monomials by canceling unwanted 
contributions. 

Some of our results are the parameterized analogues of recent work by the 
first author [21 [3] , all using a determinant summation idea that essentially 
goes back to Tutte. To make these ideas work in the parameterized setting 
is not straightforward. In particular, while the Hamiltonian path algorithm 
[3] uses determinants to cancel the contribution of unwanted labeled cycle 
covers, our /c-path algorithm from Theorem Q] uses a combinatorial argument 
to pair unwanted labeled walks of k vertices. With k = n we recover a 
0*(1.66 n ) Hamiltonian path algorithm, but using a different (and arguably 
more natural) approach. On the other hand, the parameterized packing and 
matching results of Theorems [21 [H and El and also the edge colouring result 
in Theorem [6] all use determinants. 

Our algorithm for /c-path seems to be subtle in the sense that we see no 
way of extending it to other natural combinatorial structures, even directed 
paths. In contrast, the ideas of Koutis [21] and Williams [37] work directed 
graphs, and for detecting k- vertex trees and /c-leaf spanning trees [22] in 
time 0*{2 k ). 

It seems to be difficult to achieve our results using previous techniques. 
In particular, the limitations of the group algebra framework [21j [37] were 
studied by Koutis and Williams [22]; they show that multilinear polynomials 
of degree k cannot be detected in time faster than 0*(2 k ) in their model. 



Koutis [21] has argued that the color coding method [T] and the randomized 
divide-and-conquer approach [6]) also cannot achieve running times whose 
exponent base is better than 0*(2 k ). 

2. A Projection Sieve for /c-Paths 
This section establishes Theorem [TJ 

2.1. Overview. In this section, we develop an inclusion-exclusion sieve over 
multivariate polynomials for the A;-path problem. 

The input graph is randomly partitioned into two sets V\ , V2 of roughly 
equal size. Central to our analysis is the family of labeled walks, defined 
relative to such a random partition as follows: Each occurrence on the walk 
of a vertex in G[Vi] and of an edge in G\V^\ receives a unique label. We 
call this a bijective labeling. The labels need not correspond to the order in 
which the objects are visited, and the same object can incur more than one 
label. For example, with V\ = {a, b, c}, V2 = {d, e, /}, 






Note the asymmetry between what is labeled in G\V\\ and G[V2]; indeed, 
with good probability a path of length k has roughly k/2 vertex labels in 
G[V±] but only fc/4 edge labels in G[T^]. Very roughly speaking, the running 
time of our algorithm is around 2 k ' 2+k ' A for that reason. 

With each such labeled walk we associate a monomial consisting of vari- 
ables x e for every edge e on the walk, variables y V) i for every vertex v & G[V\] 
labeled i, and variables z e> i for every edge e € G[V2] labeled i. For example, 
the monomial associated with leftmost labeled walk above is 

Xab ■ X a d ■ Xbc ' X de ■ X e f ■ y a ,i • y b ,2 ■ Vc,3 ' Zde,2 ' Z e f t l ] 

the monomial associated with the middle labeled walk contains the same xs, 
but the remaining factors are y a ,3 ■ 2/6,1 • y c ,2 • Zde,l ' z ef,2- 

If the labeled walk is a path (i.e., it has no repeated vertices), then it can 
be uniquely recovered, including the labels, from its associated monomial 
and knowledge of the source vertex. On the other hand, two different labeled 
walks that are not paths can have the same monomial, for example, 
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FIGURE 1. Representative examples of the pairing of labeled 
walks. Left: if the walk's first closed subwalk starts at a vertex 
in Vi = {a,b,c} (here, vertex c), then the labels at this vertex 
are transposed, while the walk remains unchanged. Right: if the 
walk's first closed subwalk starts at a vertex in V-z = {d, e, /} (here, 
vertex /), then the subwalk and its labels are reversed. 



In fact, we will set up a pairing such that every labeled non-path has ex- 
actly one such partner. In particular, their (identical) monomials will cancel 
each other when added in a field of characteristic 2, and only the monomials 
corresponding to labeled paths will remain. Representative examples of this 
pairing are given in Figure [TJ A good part of our exposition is devoted to 
a very careful description of this pairing; to appreciate why such caution is 
necessary, note that walks like 




are not correctly handled by our set-up and indeed will require an exception 
in the definition. 

In summary, there are four key ingredients. First, the labeled A:-walks 
that are paths will be associated with monomials that have distinct variable 
supports. Second, the monomials associated with bijectively labeled A;-walks 
that are not paths will cancel over a field of characteristic 2, establish by a 
pairing argument. Third, an inclusion-exclusion sieve will be used to cancel 
all walks that are not bijectively labeled. The sieve requires at most 3k/4 
labels (with high probability) if we are careful about the walks we consider. 
Fourth, the polynomial of A;-walks that avoid a given set of labels can be 
evaluated in time polynomial in n. 

2.2. Preliminaries on strings. From a technical perspective it will be 
convenient to view a walk in a graph as a string. To this end, let us review 
some basic terminology on strings. 



Let Abe a set whose elements we view as the symbols of an alphabet. A 
string of length £ over A is a sequence S = s\S2 ■ ■ ■ sn with Si € A for each 
i = 1, 2, . . . ,£. We say that Si is the symbol at position i of the string. The 
reverse of a string S = S1S2 ■ ■ ■ S£ is <S = s^_i ■ • • si. The concatenation of 
two strings S 1 = S1S2 • ■ ■ S£ and T = ti^ ■ • • £& is ST = S1S2 • • • S(t\t2 ■ ■ ■ ijfc- A 
string T = t\t2 ... t^ is a substring of a string 5" = S1S2 ■ • • S£ if there exists 
a j = 1, 2, ...,£ — k + 1 such that ij = Sj+j-i for all z = 1, 2, . . . , k. A 
palindrome is a string that is identical to its reverse and has length at least 
2. For Ai, A2, ■ ■ ■ , A(_ C A, we say that a string S1S2 • • • sg is an A1A2 • • • Ag- 
string if Sj £ Aj holds for every i = 1, 2, . . . , £. 

2.3. Walks in graphs. We assume that all graphs are undirected and con- 
tain neither loops nor parallel edges. For a graph G, we denote the vertex 
set of G by V = V(G), and the edge set of G by E = E(G). For convenience, 
we assume that V and E are disjoint sets. 

A k-walk in G is a string of length Ik — 1 such that 

(a) each odd position contains a vertex of G; 

(b) each even position contains an edge of G; and 

(c) for every i = 1, 2, . . . , k — 1, the edge at position 2i joins in G the 
vertices at positions 2i — 1 and 2i + 1. 

The first and last positions of a walk are the ends of the walk. 

A walk is a pa£/i if each vertex of G appears in at most one position of 
the walk. A fc-walk is closed if its ends are identical and k > 2. A walk that 
is not a path always contains at least one closed subwalk. 

Let W be a walk. A subwalk of M^ is a substring of W that is in itself 
a walk. Put otherwise, a subwalk of W is a substring with ends at odd 
positions of W. 

2.4. Representing strings as sets. Let a\02 •••o-t be a string over an 
alphabet A. It will be convenient to view a string as a set consisting of pairs 
(a,i) € A x {1,2, ...,£}, where the pair (a,i) indicates that the symbol a 
occurs at position i, that is, a^ = a. 

For a subset B C ^4 and a string a\a2- ■ ■ ae, introduce the notation 

5{aia 2 ■ ■ ■ a*} = {(a,, i) : a, e 5} C ,4 x {1, 2, . . . , £}. 
In particular, we can recover the string a\a2 ■ ■ ■ ag from the set A{a\02 ■ ■ ■ ae}. 

2.5. Admissible walks. To reduce the number of labels in the sieve, we 
will focus on a somewhat technical subset of A;- walks that we will call "ad- 
missible" walks. 

Let G be a graph with vertex set V and let s be a fixed vertex of G. 

Partition the vertex set into two disjoint sets V = V± U V2. Denote by E\ 
the set of edges of G with both ends in V\. Denote by Ei the set of edges 
of G with both ends in V2. Let k, k\,£i be nonnegative integers. 

Let us say that a fc-walk W in G is admissible if 

(a) W starts at s; 

(b) \V 1 {W}\ = k 1 ; 

(c) \E 2 {W}\ = £ 2 ; and 

(d) VF is V^-EVi-EV^-palindromeless. 
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FIGURE 2. Terminology and notation for walks introduced in 
^12. 31 12.51 The 6-vertex walk W\ is encoded by the 11-letter string 
aCdDbEeljGc. With V x = {a,b,c} we have E 2 = {H,I} and 
the "projections" Vi{Wi} = {(a, 1), (b, 5), (c, 11)} and £ 2 {Wi} = 
{(/, 8)}. The walk W\ is admissible with the following param- 
eters: starting vertex s — a, number of vertices k — 6, num- 
ber of Vi -vertices fci = 3, and number of i?2-edges £2 = 1. The 
walk W2 = elfGcGfFbAa is not admissible, because it contains the 
V^-EV^i-EV^-palindrome fGcGf. For readability, in all other exam- 
ples we denote edges by node pairs, for instance writing ef instead 
of /, with the understanding that ef is viewed as a single symbol 
for purposes of indexing W. 



Pr |Vi{P}| = % and \E 2 {P}\ = £ 2 



Here the term "palindromeless" refers to the property that a string has no 
palindrome as a substring. By V^-EVi-EV^-palindromeless we refer to the 
lack of palindromes that are also V^-EVi-EVVstrings. Observe that paths are 
palindromeless and hence V^-EVi-EWpalindromeless. 

2.6. Random projection. For a fixed ordered partition (Vi, V2), every k- 
path P in G that starts at s is admissible for some parameters k\,l<i- Con- 
versely, for fixed parameters k\,li and a fixed /c-path P that starts at s, if we 
select (X / i,V2) uniformly at random, then P is admissible with probability 
given by the following lemma. 

Lemma 7 (Admissibility). Let k\,(.2 be nonnegative integers and let P be 
a k-path in G. For (Vi, V2) selected uniformly at random, we have 

2 - k ( fci + i \A-fci-r 

^ - kx - £2) V ^2 

Proof. There are 2 fc strings of length fc over the alphabet {1,2}. The proba- 
bility in question is exactly the fraction of such strings that have exactly k\ 
1-positions and exactly £2 22-substrings. There are exactly k\ + 1 positions 
where to interleave the k\ Is with substrings of 2s. Each such substring of 
length j contributes exactly j — 1 22-substrings. The total number of 2s is 
k — ki, so there must be k — k\ — £2 substrings of 2s. The positions where the 
substrings interleave the Is are allocated by the first binomial coefficient. It 
remains to allocate the lengths of the strings. The total length is k — k%, 
and each of the k — k\ — £2 strings must have length at least 1 . Thus there 
are k — ki — {k — k\— £2) = £2 free 2s to allocate to k — k\ — £2 distinct bins. 
The second binomial coefficient carries out this allocation. □ 

In particular, a fixed fc-path starting at s is admissible with positive prob- 
ability if and only if either k\ = k and £2 = or k\ < k and k\ + £2 < k — \ < 

2k t +£ 2 . 
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Let us now derive an asymptotic approximation for the probability in 
Lemma [JJ We employ the following variant of Stirling's formula due to 
Robbins [32J. For all j = 1,2, . . . it holds that 




(1) j! = v / 2vri( ) e<Ej where -,„.,., < £j < 



1 1 

12j + 1 3 Uj 



Let us abbreviate 



t) = (^Y( l - b a 



From Stirling's formula ([T]) it follows that (?) = ©*((^)) holds uniformly for 
all < b < a < n. We can thus approximate the probability in Lemma [7J 
uniformly for < £ 2 , k\ < k such that k\ + £2 < k — 1 < 2k\ + £2, with 



Pr( \Vi{P}\ -=jfci and |£ 2 {.P}| 
(2) 



e * (2 -,/ fei \A-fci 



k-h-£ 2 

2.7. Labeled admissible walks. The following labeling scheme for admis- 
sible walks serves two purposes. First, labeling enables us to "decouple" the 
sieve from the graphical domain (that is, vertices and edges) into a set of 
abstract labels whose number depends only on the parameters k\,ti and 
not on the size of the graph. Second, the labeling facilitates cancellation of 
non-paths in the sieve. 

Let K\ be a set of k\ labels. Let L2 be a set of £2 labels. For example, 
let K x = {1,2,..., A*} andL 2 = {l,2,...,^ 2 }. 

Let W be an admissible walk. Let k\ : Vi{W} — ► K\ and A 2 : -E 2 {iy} — \ 
L 2 be arbitrary functions. The three-tuple (W, «i, A 2 ) is a labeled admissible 
walk. Intuitively, each position in W that contains a vertex in V\ gets 
assigned a label in K\ by K\. Similarly, each position in W that contains an 
edge in E2 gets assigned a label in L 2 by A 2 . Let us say that the labeling is 
bijective if both k\ and A 2 are bijections. 

Example. Consider two labelings of the same walk W, 




The walk is admissible with parameters s = c, k = 6, k\ = 3, and £2 = 
2. Both labelings associate a label with each position of W that con- 
tains a symbol from V\ = {a,b,c} or E 2 = {de,ef}. We have T4{W} = 
{(6, 3), (6, 11), (c, 1)}, that is, there are three occurrences of symbols from V\ 
in W; in particular the symbol b occurs at the 3rd and the 11th position. 
Similarly, we have E^jVF} = {(e/, 6), (de, 8)}. The labeling on the left is 

Ai(6,3) = 2, Ai(Ml) = l, Ai(6,l) = 3, «2(e/,6) = 2, « 2 (de,8) = l. 

We observe that this labeling is bijective since Ai and k 2 are bijections. 



12 

The labeling on the right is 

Ai(6,3) = 3, Ai(6,ll) =Ai(c,l) = l, K 2 (ef,6) = K 2 (de,8) = 2, 
and not bijective. In fact, Ai avoids the label 2, and k 2 avoids the label 1. 

2.8. Fingerprinting and identifiability. We associate with each labeled 
admissible walk an algebraic object (or "fingerprint") that we use to repre- 
sent the labeled admissible walk in sieving. Here it is important to observe 
that while we are careful to design the fingerprint so that each labeled path 
has a unique fingerprint, the fingerprints of labeled non-paths are by de- 
sign not unique — we will explicitly take advantage of this property when 
canceling labeled non-paths in N2.101 

The sieve operates over a multivariate polynomial ring with the coeffi- 
cient field ¥ 2 b (the finite field of order 2 ) and the following indeterminates. 
Introduce one indeterminate x e for each edge e G E. Introduce one indeter- 
minate y v j for each pair (v,i) G V\ X K\. Introduce one indeterminate z e ^ 
for each pair (e,i) G E 2 x L 2 . 

Let (W,ki,X 2 ) be a labeled admissible walk. Associate with (W, n±,X 2 ) 
the monomial (fingerprint) 

m(w,Ki,x 2 )= n Xe n vv,m(!>A n ^m- 

(e,j)eE{W} (v,i)eVi{W} (e,i)£E 2 {W} 

The following lemma is immediate. 

Lemma 8 (Identifiability). The monomial m(W, k,\ , X 2 ) of a labeled admis- 
sible walk (W,ki,X 2 ) uniquely determines the edges and their multiplicities 
of occurrence inW. In particular, any path is uniquely identified. Further- 
more, if W is a path and K\,X 2 are bijections, then m(W,Ki,X 2 ) uniquely 
identifies (W, K\,X 2 ). 

Example. We presented some example monomials already in £12.11 We can 
also consider a bijectively labeled walk that repeats an edge in E 2 : 




x le ■ x bf ■ xlf ■ Ub,l ■ Vb,2 ' Z e f,l ■ Z e f >2 , 



d e '-' / 

and a non-bijectively labeled walk, 



H3 




2 H^p 2 i 

e 



/ 



Xbc ■ x bd ■ x bf ■ x de ■ x ef ■ y b)1 ■ y 6i3 • y Ctl ■ z de ^ 2 ■ z e f >2 . 



In all these examples observe that if the walk is a path, we can reconstruct 
it from the x-variables and knowledge of the start vertex. Because a path has 
neither repeated vertices nor edges, the y- and z-variables in the monomial 
enable us to reconstruct the labeling. 
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2.9. Sieving for bijective labelings. Let us denote by £ the set of all 
labeled admissible walks. For I\ C K\ and J 2 C L 2 , denote by £[Ii, J 2 ] the 
set of all labeled admissible walks that avoid the labels in I\ and J 2 . Let us 
denote by 33 the set of all bijectively labeled admissible walks. 

By the principle of inclusion-exclusion, we have 

(3) Yl rn{W,K 1 ,X 2 )= £ £ (-1)W+I J »I £ m(W,K 1 ,X 2 ). 

{W,mM)e& ZiC^jJaCLa (W,«i > A 2 )s£[Zi,J 2 ] 

2.10. Bijectively labeled non-path fingerprints cancel. Let us parti- 
tion 33 into 33 = 7 U 31, where CP consists of bijectively labeled admissible 
paths, and 31 consists of bijectively labeled admissible non-paths. Accord- 
ingly, the left-hand side of ([3]) splits into 

Y m(W,K 1 ,X 2 )= J2 rn(W,Ki,\ 2 )+ Yl m(W, m, X 2 ) . 

(W,«i,A 2 )63 (W,Ki,A 2 )e:P (W,«i,A 2 )eK 

We show that the rightmost sum vanishes. To this end, let us first re- 
call that an involution is a permutation that is its own inverse. We claim 
that it suffices to construct a fixed-point-free involution : 3i — >• 3i with 
m(W, «i, A2) = m,(</)(iy, «i, A2)) for all (W, «i, X 2 ) G 31, Indeed, introduce an 
arbitrary total order to 31 and observe that in characteristic 2, we have 

Y m(W,Ki,\ 2 )= Y rn(W,Ki,\ 2 ) + m(<l>(W,K 1 ,\2)) 

(W,re 1 ,A 2 )eK (W,rei,A 2 )6S 

(W,k 1 ,A 2 )<^(W,«i,A 2 ) 

2m(W,Ki,X 2 ) = 0. 



E 



(TV,ki,A 2 )6S 

(iy,Ki,A 2 )<0(w,Ki,A 2 ) 

To construct a fixed-point-free involution : 3J — > 31 with m(W, m, X 2 ) = 
m(4>(W, «i, A2)) for all (W, «i, A2) € 3?, we observe that every walk W that 
is not a path contains at least one closed subwalk. In particular, W contains 
a first closed subwalk, that is, the closed subwalk C with the property that 
C is the unique closed subwalk in the prefix SC of W = SCT. 

We denote the first closed subwalk of W by C(W) and by c{W) the first 
(and hence also the last) vertex of C(W). 

Let us partition 31 into two disjoint sets, 3l\ and 3i 2 , where 

3i 1 = {(W,K 1 ,X 2 ) &3l:c(W)eV 1 }, 

3i 2 = {(W,K 1 ,X 2 ) G 31 : c(W) G V 2 } . 

We proceed to construct the pairing 4> on these two sets. See Figure Q] for 
examples. 

2.11. The pairing on 3l\ — label transposition. Select an arbitrary 
(W, Ki,X 2 ) G 3l\. Let j and £ be the positions of W that contain the symbol 
c(W) and constitute the ends of C(W). For brevity, let us write c for c(W). 
Because c G V\, we have (c,j), (c,t) G Vi{W}. Define k[ to be identical to 
K\ except that 

«i (c, i) = «i (c, £) , «i (c, £) = K\ (c, j) . 

Observe that K±(c,j) ^ K\(c,£) because (W,Ki,X 2 ) is bijectively labeled. 
Thus, k[ 7^ K\ and (W,«/ 1 ,A2) G 3?i. Furthermore, we have m{W,K,\,X 2 ) = 



14 

m(W, k[, A2). Thus, we can set (j>(W, K\, \ 2 ) = (W, k[, X 2 ) to obtain the de- 
sired fixed-point-free involution on 3?i. Indeed, <^(W, «i, A2) = (W, k[, A2) 7^ 
(W, «i, A 2 ) and 2 (^, «i, A 2 ) = (W, «i, A 2 ). 

2.12. The pairing on 3?2 labeled reversal of first closed subwalk. 

Select an arbitrary (W, m, A2) € 3?2- Let C = C(VF) and let S 1 , T be strings 
such that 

W = SCT . 

Let us define the string W by reversing C in W, that is, 

W' = stJT . 

We observe that the strings C and C have identical ends because C is a 
closed walk, implying that W' is a walk in G. We also observe that W is 
admissible. Indeed, because c{W') = c{W) G V2, any V^Vi-EV^-palindrome 

ueveu in W' can either (a) occur as a subwalk of C in W', or (b) have at 
most one position in W' common with C . But in both cases we have that 
ueveu occurs in W, which is a contradiction since W is admissible. Thus, 
W is admissible. 

We observe that W = W if and only if C is a palindrome. Furthermore, 
if we reverse C(W) = C in W', we obtain back W. That is, W" = W. 

In terms of string positions, we can characterize the reversal W i-> W' 
using the following permutation of positions. Let j and £ be the positions of 
W that constitute the ends of C. Define the permutation p : {1, 2, . . . , k} — >■ 
{l,2,...,/c}by 

I i if z < 7 or i > £; and 

l£-i+i i£j<i<e. 

Let us denote the symbol at the ith position of W by u>i. The reversal 
W i-> W' can now be characterized by observing that w'u\ = w% holds for 
each i = 1, 2, . . . , k. 

We now introduce a labeling k[, A 2 of W' using the labeling m, X 2 of W. 
In particular, let us label W so that each position of W is labeled using the 
label of the /^-corresponding position in W, if any. In precise terms, using the 
labeling m : Villi 7 } — >■ K\, define the labeling k[ : Vi{W'} -> K\ for each 
(«^ (i) ,p(i)) G MW} by setting k'^w'^, p(i)) = m(wi,i). Similarly, using 
A 2 : E 2 {W} -> L 2 , define A' 2 : E 2 {W'} -> L 2 by setting A' 2 (u£ (i) ,p(i)) = 
A 2 («;»,«) for each (u^ (i) ,p(z)) e £ 2 {W}. 

Now set (f)(W,m, A2) = (W',Ki,A 2 ) and observe that 0(W, «i, A2) € #2, 
2 (W,ki,A 2 ) = (W,ki,X 2 ), and m(W, «i,A 2 ) = m(0(lV, m, A 2 )). 

What is not immediate, however, is that 4>(W, m, A2) 7^ (W, «i, A 2 ). There 
are two cases to consider, depending on C. 

In the first case, C is not a palindrome, that is, C 7= C. Thus, W' 7= W 
and hence (W, n[, A 2 ) / (H 7 , «i, A 2 ). 

In the second case, C is a palindrome. Since C is a closed walk, the 
string C has odd length at least 3. In particular, the length 3 (that is, a 
palindrome of the form ueu with u S V 2 and e £ E) cannot occur because 
G has no loop edges. For palindromes of length 5, the only possibility is 
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that C is a V^ii^V^-E^V^-palindrome. Indeed, C can neither be a V\EV\EV\- 
palindrome nor a yi-EV 2 £'Vi-palindrome because c(W) € V2. Furthermore, 
C cannot be a V^-EVi-EV^-paiindrome because such palindromes by def- 
inition do not occur in the admissible W . Thus, for length 5 the only 
possibility is a V^E'V^-EV'rpalindrome, that is, a V^-E^V^E^V^-palindrome. 
Such a palindrome contains two occurrences of an edge in E2 that are in 
p-corresponding positions. These occurrences get different labels under A2 
and A 2 . Thus, (W, k[, A' 2 ) 7^ (W, «i, A2). Finally, we observe that C cannot 
have length more than 5, because a palindrome of length 7 or more must 
include a palindrome of length 5, which would contradict the assumption 
that C is the first closed subwalk in W. 

2.13. The algorithm. First, we recall the following result: 

Lemma 9 (DeMiilo-Lipton-Schwartz-Zippel [91 [33]). Let p(xi, x 2 , • • • , x n ) 
be a nonzero polynomial of total degree at most d over the finite field ¥ q . 
Then, for 0,1,0,2, ■■■,a n € ¥ q selected independently and uniformly at ran- 
dom, 

Pr(p(ai,a 2 ,...,a n ) / 0) > 1 . 

q 

Let us assume the parameters k, fci,-£ 2 have been fixed so that k\ + £2 < 
k — 1 < 2k\ + £2- (We will set the precise values of fci,^ 2 in what follows.) 
To decide the existence of a /c-path starting at s, we repeat the following 
randomized procedure. 

First, the procedure selects an ordered partition (Vi,V2) uniformly at 
random among all the 2 n such partitions. Lemma [7] implies that a fixed 
A;-path P that starts at s is admissible with positive probability. 

Next, the procedure makes use of Lemma [9] to witness a nonzero eval- 
uation of a multivariate generating function for labeled admissible /c-paths 
starting at s. In particular, from ([3]) and ^2.10l we have that 

(4) Yl rn(W, Kl ,X 2 )= £ £ (-lfil+^l ^ m(W,«i,A 2 ). 

(W / ,Ki,A2)eJ' ZiCXiJaCLa (VK,ki,A 2 )g£[/i,J 2 ] 

The left-hand side of ([U is a multivariate polynomial of degree at most 
k — l + fci+^2- It follows from Lemma[8]that the polynomial is not identically 
zero if and only if G has an admissible /c-path starting at s. 

It remains to evaluate the right-hand side of (JU) for a random assignment 
of values in F 2 t to the indeterminates. To this end, the procedure iterates 
over each I\ C K\ and J 2 C L 2 and employs dynamic programming to 
evaluate the rightmost sum in (j3|). 

Without loss of generality we can assume k > 3. For parameters k, k%,£2 
and a string T = £1^2^3^4*5 over the alphabet V U E, our objective is to 
compute 

(5) M(k,k 1 ,£ 2 ,T)= J2 m(W,K 1 ,X 2 ). 

T is a suffix of W 

In particular, taking the sum over all T, we obtain the rightmost sum in (j3|). 

The recursion for ([5]) is as follows. For a logical proposition P, let us 

define [P] to be 1 if P is true and otherwise. For k > 3, we observe by 
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induction on k that 

M(k,k 1 ,e 2 ,t 1 t 2 t 3 t 4 t 5 ) = 

= [^1^2*3*4*5 is not a V^-EVi-EV^-palindrome] 

(6) x J2 (b i Vi] + [« € Vi] ^ y^-) ([e £ £ 2 ] + [e e £ 2 ] £ z ej ) 

ees jei^i\/i jei 2 \J 2 

e = {Mi} 

x Af(& -Ui-[»£ Vl],£ 2 - [e € £ 2 ],wii£ 2 £ 3 ) • 

To set up the base cases for the recursion, we observe that M(k, ki,£ 2 ,T) 
can be computed for all < ki,£2 < k = 3 and all T = ^1*2*3*4*5 in time 
polynomial in n. Furthermore, M(k,k\,£ 2 ,T) = whenever k\ > k or 
£ 2 > k or h\ < or £ 2 < 0. 

Consequently, for any given assignment of values in ¥ 2 b to the indetermi- 
nates x e , y v> £, and z e i, the procedure evaluates the right-hand side of dU 
via © in 0(2 kl+i2 k 3 n 4 ) arithmetic operations over ¥ 2 b. 

Let us now complete the algorithm by optimizing the parameters for run- 
ning time and 0(1) probability of success. Denoting the probability that a 
fc-path P starting at s is admissible by P(k, k\,£ 2 ), we have that in r repeti- 
tions of the procedure at least one repetition finds P admissible with prob- 
ability 1 - (1 - P(k, hJ^Y > l-e~ p ^ k ^>. Settingr = \1/P(k,k u £ 2 )] 
and b = |~log 2 6k] , it follows from Lemma [9] that any fixed /c-path start- 
ing at s in G is witnessed with probability at least (1 — e _1 )/2 in time 
0(2 kl+£2 k 5 n 4 /P(k,k 1 ,£)). Setting ki = [-fik\, £ 2 = [-f 2 k\, and employing 
([2]) to approximate P(k, k±,£ 2 ), we obtain 0*(1.6569 ) time for 71 = 0.5 and 
72 = 0.207107. 

3. A Determinant Sieve for ^-Dimensional p-Packings 
This section establishes Theorem [SJ 

3.1. Prepackings and Edmonds's symbolic determinant. Let us say 

that a subset A C 3" is a j-prepacking if \A\ = j and the sets in A are 
pairwise disjoint when projected to U\ U U 2 . 

Observe that each A £ A in a j-prepacking identifies both a unique 
ui(A) £ A n XJ\ and a unique u 2 {A) £ At~)U 2 . 

For a bijection a : U\ — > U 2 , let us say that a j-prepacking A is compatible 
with a if for all A £ A it holds that a(u\{A)) = u 2 (A). Note that each j- 
prepacking is compatible with at least one a. 

Edmonds |10j made the algorithmically seminal observation that the de- 
terminant of a symbolic r x r matrix E = (e MljU2 ) Ulg (7 liM2g (7 2 is a signed 
generating function over partitions of U\ U U 2 into 2-subsets with exactly 
one element from U\ and exactly one element from U 2 . Indeed, identifying 
each such partition with a bijection a : U\ — )• U 2 , we have 

(7) detE= ^2 s g n r(o-) II e«i,a(«i)> 

<T:t/i->(72 «ie£/i 

<T bijective 

where the sign sgn T (<j) is the sign of the permutation gt for an arbitrary 
fixed bijection r : U 2 — > U\. 
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Our strategy is to leverage Edmonds's observation from the dimensions 
U\ and U2 into q dimensions U±, U2, ■ ■ ■ , U q with sieving. In particular, Ed- 
monds's observation forces the packing constraint in the first two dimensions, 
which allows us to restrict the sieve to the remaining q — 2 dimensions. 

3.2. Fingerprinting and identifiability. Consider a j-prepacking A C 3~. 
The domain of the prepacking is the set 

(8) d(A) = {{u,A) : u G A G A} C (C/ 3 U U 4 U • • • U U q ) x J. 

Observe that \d(A)\ = j(q - 2). 

Let L be a set of p(q — 2) labels. A labeling of .A is a pair (<r, A), where 
a : U\ — > C/2 is a bijection compatible with A and A : <i(yl) — )• L is an 
arbitrary mapping. The labeling is bijective if A is a bijection. We say that 
a triple (A, a, A) is a labeled j-prepacking. 

The sieve operates over a multivariate polynomial ring with the coefficient 
field F 2 f> and the following indeterminates. Introduce the indeterminate 
w for tracking the weight j of a j-prepacking. Associate with each A S 
3" an indeterminate xa- Associate with each pair (7/1,^2) £ U\ x C/2 an 
indeterminate y Ul ,u 2 - Associate with each pair (u,£) G (L^UL^U- • -UUg) xL 
an indeterminate z Ui £. 

The signed monomial of a labeled j'-prepacking (A, a, A) is 

(9) m(A,a,X) = sgn T (a)w 1 Y[x A Y[ y Ul ,a( Ul ) \\ z u,\{u,A) ■ 

AeA utGU (u,A)ed(A) 

Lemma 10 (Identifiability). The monomial m(A, a, A) uniquely determines 
both A and a. Furthermore, if A is a p-packing and A is bijective, then 
m(A,o~,X) uniquely determines A. 

3.3. Sieving for bijective labelings. Denote by L the set of all labeled p- 
prepackings. For J C L, denote by £[J] the subset of labeled p-prepackings 
whose labeling avoids each label in J. Denote by ¥> the set of all bijectively 
labeled p-prepackings. 

By the principle of inclusion-exclusion, 

(10) J2 m(A,a,A) = ^(-l)l J l Yl rn(A,a,X). 

(A,<r,A)eS JCL (A,a,\)eL[J] 

3.4. Fingerprints of bijectively labeled non-p-packings cancel. Par- 
tition S into "B = TWJl, where 7 is the set of bijectively labeled p-packings, 
and R is the set of bijectively labeled p-prepackings that are not packings. 
Accordingly, the left-hand side of (jlOp splits into 

y m(A,a, A) = \, rn(A,a,X)+ \, m(A,a,X). 

{A,a,\)e-B (A,cr,\)&> (A,(T,\)&k 

We show that the rightmost sum vanishes in characteristic 2. To this end, 
it suffices to construct a fixed-point-free involution <p : 31 — > 31 such that 
m((f)(A, a, A)) = m(A, a, A) holds for all (A, a, A) G 31. Consider an arbitrary 
(A, a, X) G 31. Since A is a p-prepacking but not a packing, there is a 
minimum (with respect to e.g. lexicographic order) three-tuple (uq, A±,A2) G 
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(U s U Ui U • • • U U q ) x A x A such that u e A 1 n A 2 and Ai / A 2 . Define a 
labeling A' : d(A) — )• L of .A by setting, for each (u, A) € d(A), 



(11) A'(u,A) 



A(u, A) if «^«o or A^ {Ai,^ 2 }; 
X(uq, A2) if n = no and A = A\\ and 
X(uq, -Ai) if u = uq and A = A 2 . 



Note that A' is bijective and that A' 7^ A. From (fl8|) and (fTT|) it follows 
that m(yi, 0", A') = m(A, a, A) holds for all (A, a, A) € 3?. We can now 
set (p(A, a, A) = (A, a, A') and observe that 4>(A, a, A) € 3?, <^>(.A, ex, A) / 
(.A, <7, A), and 4> 2 (A,a,\) = (.A,^, A) for all (A, a, A) € 3?. Thus, is a 
fixed-point-free involution on 31. 

3.5. The algorithm. From (fl"0j) and £ )3.4l we have 

(12) 2 m(A,a,A) = ^(-l)l J l £ m(A«7,A). 

(A,a,A)ey JCL (A,A)e£[J] 

From Q we observe that the left-hand side of (|12|) is a multivariate 
polynomial of degree at most pq + r. It follows from Lemma [TU] that the 
polynomial is not identically zero if and only if 3~ contains a p-packing. It 
remains to evaluate (|12|) for an assignment of values to the indeterminates. 

Let JCLbe fixed. Introduce an r x r matrix E( J) as follows. Index the 
rows with U\ and the columns with XJi. Define the entry at row u\ £ U\, 
column u 2 € U2 by 

e Ul ,u 2 (J) =y Ul ,u 2 { 1 + w ^] ^ [[ E w)- 

Denote by L j [J] the subset of labeled j-prepackings whose labeling avoids 
each label in J. From ((TJ), {(HJ) , and @ it immediately follows that we have 

r 

(13) det E(J) = J2 Yl rn(A,a,X). 

j=0 {A,a,\)&Cj[J] 

Consequently, for any given assignment of values in F 2 b to the indetermi- 
nates xa, Vu 1 ,u 2 i an d z u,e, we can evaluate the left-hand side of (fl3j) as a 
polynomial in the indeterminate w via (|13j) . Taking the sum over all J C L 
and extracting the coefficient of the monomial w p , we obtain an evaluation 
of the right-hand side of (|12p in total 0(2 p ( q ~ 2 '\3^\pq 2 r 4 ) arithmetic opera- 
tions in F 2 6. Taking b = |~log 2 2(pq + r)~|, Lemma [9] implies that we witness 
a p-packing in 3" (as a nonzero evaluation of the left-hand side of (|12p ) with 
probability at least 1/2 in time O* (2 p ( q ~ 2 >) for polynomial size families 3". 



4. A Projection-Determinant Sieve for ^-Packings of <?-Sets 
This section establishes Theorems [2] and HI 
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4.1. Tutte's observation. Let us recall that an involution is a permutation 
that is identical to its inverse. In particular, the cycle decomposition of an 
involution consists of fixed points and transpositions (cycles of length 2). It 
follows that the involutions on a set I are in a one-to-one correspondence 
with the set partitions of I into sets of cardinality 1 and 2. The following 
lemma is essentially due to Tutte 



Lemma 11 (Tutte's Determinant-Partition Lemma). Let T be an m x m 
matrix with entries in a multivariate polynomial ring over a field of charac- 
teristic 2. Index the rows and columns of T by the elements of a set I and 
suppose that T is symmetric so that 



t*3 



z2k s {i},k tfi — j; 



holds for all i,j € I. Then, 



( 14 ) detT = J2 II E S W- 

l-.I-^I iel k 

i involution i<i{i) 

Proof. Denote the set of all permutations of I by Pj. Observe that a per- 
mutation a 6 Pi is not an involution if and only if the cycle decomposition 
of a contains a cycle of length at least 3. Suppose that v G Pi is a per- 
mutation that is not an involution. Introduce an arbitrary total order on 
/, and order the cycles of length at least 3 in v based on the least point in 
/ moved by each such cycle. Denote by v' the permutation obtained from 
v by inverting the first cycle of length at least 3 in v. Clearly, v' ^ v and 
(i/y = v. Now observe that because T is a symmetric matrix, for a cyclic 
permutation [i\ ii ■ ■ ■ ij) and its inverse (ij i,_i • • • i\), we have 



Hl,i2^2)*3 ' ' ' ij—ltijijiH ij4j—iij—i4j—2 ' ' ' *2)*iHi, 



It follows that 






Partition Pi into Pi = Qi U Ri, where Qi consists of the involutions 
and Ri consists of the non-involutions. Introduce an arbitrary total order 
on Ri. Because the determinant of T is equal to the permanent of T in 
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characteristic 2, we have 

det T = J2YI V(0 
ctGPj iei 

= y n*v(o+ iz n **."(*) 

tGQ/ ie/ 1/gJ?; iei 

= y n **.*(*) + ^ (nvw+nvw) 
= y n **.*(*) + ^ 2 n *<."W 

ieQi iei veRi iei 

v<v' 

= y n*vw 

Thus, splitting the product over fixed and moved points of i, and using 
the symmetry of T, we have 

det T = Y II *».<•(») II *»>'(») 

teQj 16/ 16/ 

j=t(i) i^i(i) 

y n E4}.* n (S 5 {vw}.fc 

tGQ/ 16/ fc iei k 

= y n E s «,fc n (5Z 5 {<},fc +2 S s {vw}.* 5 {v(0}.*') 

i.eQi iei k iei ^ k k<k' ' 

i=i(i) i<i(i) 

= y n E a {i},* n Y s iiAi»,k 

ieQi iei k iei k 

i=i(i) i<c(i) 

= Y II Y4iAi)},k- 

teQi iei k 

i<i{i) 

The claim follows. D 

Our strategy is to leverage Tutte's observation with random projection 
and sieving. In particular, we witness a p-packing by randomly projecting it 
to a set U\ QU where Tutte's observation forces the packing constraint with 
positive probability, which allows us to restrict sieving to the complementary 
projection into U2 = U \U\. 

4.2. Admissible packings and prepackings. Let 3" be a set of g-subsets 
of an n-element universe U. Partition U into two disjoint sets U = U± U Ui 
with |C/i| = n\ and IC/2I = ^2 = n — n\. We say that such an ordered 
partition (U\,Uz) of U is an (m, 722)- partition. 

A subset A C 3~ is a p-packing if |A| = p and the sets in .A are pairwise 
disjoint. We say that A is admissible if every set A € A satisfies |.An£/i| < 2. 
We say that A is a (po, p±, p2)-prepacking if 

(a) \A\ =Po+p\ +P2; 
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(b) \{A € A : \A n Ui\ = j}\ = pj for each j = 0, 1, 2; and 

(c) the sets in A are pairwise disjoint when projected to U\. 

Note that a prepacking is by definition admissible. Also, every admissible 
p-packing is a (po,.Pi>P2)-prepacking for some parameters po + p\ + P2 = P- 
In this case we say that the p-packing is a (po, pi, p2)-packing. 

Let us say that a (po;Pij.P2)-prepacking A is compatible with an involution 
l : U\ — > U\ if for every A € A it holds that 

(a) l fixes the point in A D C7i if \A D E/i| = 1; and 

(b) i transposes the two points in A n £/i if | A n U\ \ = 2. 

Note that every prepacking is compatible with at least one involution. 

4.3. Random projection. We analyze the probability that a given p- 
packing projects under a random (Ui, U2) into a (pojPi>P2)-packing. 

Lemma 12 (Admissibility). Let A be ap-packing. For an (711,712) -partition 
(Ui,U2) ofU selected uniformly at random, we have 

Yt(A is a (po,pi,p2) -packing) = 

( 15 ) = ( P \( Pl+P2 )( q ) Pl ( q Y 2 ( n ~ PQ )( n 

\Pi +P2J \ P2 J V 1 / V 2 / V n i ~Pi- 2 P2j \n\ 

Proof. Among the p pairwise disjoint sets in A, there are (, ? ) ways to 
select the sets that intersect U\ in 1 or 2 points, and (P 1+P2 \ ways to select 

among these the P2 sets that intersect in 2 points. There are (?) 1 (f\ 
possible intersection patterns with U% in these selected p\ + P2 sets. There 
are ( n ™~ v \ ) ways to select the remaining n\ —p\ — 2p2 points of U\ outside 
the pq points of A. □ 

Remark. A nonzero probability is allocated if and only if pq < n, po + 
P1+P2 = P, Pi + 2^2 < 77i, and n\ - p\ - 2p 2 < n - pq. 

Using techniques similar to £ 12.61 let us derive an asymptotic approxima- 
tion for (|15p . One verifies by direct calculation that for 5 = {p\ +2p2)/(pq), 

n — pq \ / n \ / pq 



x 5n — p\ — 2p2/ \8nj \p\ + 2^2/ 

For n\ = [5n\ thus, uniformly for all < Po,Pi,P2 < P with P0+P1+P2 = P, 
Pt(A is a (po,pi,p2)-packmg) = 

^ Q J( P \/Pl+P2\(q\(q\/ PQ V 

\\pi+P2/\ P2 /WWV + w 

4.4. Fingerprinting and identifiability. Let A C 3" be a {po,Pi,P2)- 
prepacking. The domain of the prepacking is the set 

(17) d(A) = {(u,A) :u£AeA} C U 2 x J. 

Observe that |d(.A)| = qpo + (q — l)pi + (q — 2)p2- 

Let L be a set of qpo + (q — l)pi + (q — 2)p2 labels. A labeling of A is 
a pair (t, A), where l : C/i — >■ U\ is an involution compatible with A and 
A : d(A) — > L is an arbitrary mapping. The labeling is bijective if A is a 
bijection. We say that a triple (A,t,X) is a labeled (po, pi,p2)-prepacking. 
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The sieve operates over a multivariate polynomial ring with the coefficient 
field ¥ 2 b and the following indeterminates. Introduce the indeterminates 
wo, wi, and w^ for tracking the parameters Po,Pi,P2 of A. Associate with 
each A £ 3" an indeterminate xa- Associate with each set K C U\ of size 
1 < | if | < 2 an indeterminate yx- Associate with each pair (u,£) G^xL 
an indeterminate z u £. 

The monomial of a labeled (po,pi,p2)-prep&ckmg (A, i, A) is 

, 1R . m(Aa)=^ p °^ 1 ^ 2 n 4 n 4,«} n <a ( m>- 

^ i °' ) Ae.A iel/i («,A)ed(yi) 

i<t(i) 

Lemma 13 (Identifiability) . T/ie monomial m(A, l, A) uniquely determines 
both A and t. Furthermore, if A is a p-packing and A is bijective, then 
m(A,L,\) uniquely determines A. 

4.5. Sieving for bijective labelings. Denote by £p ,p 1 p 2 the set of all 

labeled (po^Pij^-pre-packings. For J C L, denote by £p , Pl , P2 [</] the sub- 
set of labeled (po,PiiP2)-prepackings whose labeling avoids each label in J. 
Denote by 23 P0)PliP2 the set of all bijectively labeled (po,Pi,P2)-prepackings. 
By the principle of inclusion-exclusion, 

(19) J2 m(yUA) = ^(-l)l J l E m(A,L,X). 

(A,L,X)€'Bp ,p 1 ,p 2 JCL {■A-,i-,y)££'p ,p 1 ,P2[J} 

4.6. Fingerprints of bijectively labeled non-(po,Pi, J>2)-packings can- 
cel. Let 23 P0iPl)P2 be the set of bijectively labeled (po,Pi>£>2)-prepackings. 
partition &po,pi,p2 mxo -Dpo,pi,p2 — "^po,pi,P2 ^ poiPi>P2> wnere •'po>piiP2 is tne 
set of bijectively labeled (po,Pi, J>2)-packings, and Rp , Pl ,p 2 is the set of bi- 
jectively labeled (po 5 Pi)P2)-pre-packings that are not (poiPi>£>2)-packings. 
Accordingly, we have 

y m(A, o, X) = 2_j m(A,L,X)+ \, m(A,L,X). 

(^l,i,A)G-Dpg 1 p-^ 7 p2 (■^-i'')'^J^ i ^P0,Pi 1 P2 \'^) t 5^J^="^P0)Pl jP2 

By a pairing argument essentially identical to the one given in £ )3.4l the 
rightmost sum vanishes in characteristic 2. 

4.7. The algorithm. Let us assume that the parameters < po,Pi,P2 < P 
and ni,ri2 have been fixed so that any given p-packing is a (po,Pi,P2)- 
packing with positive probability. (We will set the precise values in what 
follows.) The algorithm repeats the following randomized procedure. 

First, the procedure selects an ordered (ni,ri2)-partition (f/i,^) uni- 
formly at random among all the ( " ) such partitions. 

Next, the procedure evaluates the following generating function for a ran- 
dom assignment of values to the indeterminates. From (|19p and §4.61 we 
have 

(20) Yl rn(A,L,X) = ^2(-l)^ £ m(A,L,X). 

(A,L,X)€7 POtPltP2 JQL (A,L,\)€Cp 0i p liP2 [J] 



The left-hand side of (|20p is a multivariate polynomial of degree at most 
2ni + (2q + 4)p + (2<? + 3)pi + (2q + 2)p 2 - It follows from Lemma [H that 
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the polynomial is not identically zero if and only if 3" contains a (po,pi,P2)~ 
packing. 

It remains to evaluate the right-hand side of (|20p . Let J C L be fixed. 
The procedure relies on the following observation that a sum over labeled 
prepackings factors into a product of two independent expressions. The first 
expression, S(J), generates the sets that do not intersect U\ with a simple 
product. The second expression, detT(J), generates the sets that intersect 
U\ with Lemma [TTJ 

In precise terms, let 

(21) s(j) = n (i + w ix\ n e z ki) ■ 

,401/1=0 

Define the symmetric n\ x n\ matrix T{J) as follows. Index the rows and 
columns by elements of U\. For u\ G U\, define the diagonal entries by 



(22) 



tu 1 ,u 1 (J) = y 2 {ui} U + wj E x a n E Z U2,A 



A& u 2 £AnU 2 t£L\J 

Ar\Ui={ui} 



For u\ , V\ G C/i with u\ ^ v \ , define the off-diagonal entries by 

(23) t uuvl (J) =y {uitVl} (l + w 2 E x A Y[ E z U2 A 



AnC/i={ui,ui} 



Let us now observe that 



( 24 ) E E E "^ *> A ) = s ( J ) det T ( J ) • 

0<P0<|31 0<pi+2p 2 <ni (yt,t,A)G-Cp , P1 ,p 2 [J] 

To this end, first recall (|18|) and the notion of compatibility between a 
prepacking and an involution ( §4.2p . Next, expand ()22[) and (j23|) to sums 
of monomials and apply Lemma [11] to conclude that detT(J) is exactly the 
left-hand side of (|24p restricted to po = 0. Finally, expand ()2ip to conclude 
that HMD holds. 

Consequently, for any given assignment of values in F 2 b to the indetermi- 
nates xa, Vk, and z u> i, the procedure evaluates the left-hand side of (|24"|) as 
a polynomial in the indeterminates wq, W\,W2 using a total of 0(|3"| 6 g|L|nf ) 
arithmetic operations in F 2 &. From such an evaluation we can recover the 
coefficient of the monomial w^w^w^ 2 - This coefficient corresponds to an 
evaluation of the inner sum in the right-hand side of (|20p . Taking the sum 
over J C L (and multiplying by w^u^wff), we obtain an evaluation of the 
right-hand side of (|20p . 

Denoting the probability that a p-packing A is a (pQ,pi,p2)-p&cking with 
P(n,ni,po,pi,p2), and taking r = \l/P(n,ni,po,pi,p2)~\ repetitions of the 
procedure with b = |"log 2 16n], Lemma [9] implies that at least one repetition 
of the procedure witnesses any fixed p-packing A (as a nonzero evaluation 
of ([20])) with probability at least (1 — e _1 )/2 in time 

O(2 9P0+ ( q -V pi+ ( q - 2 ') p2 \$\ 6 pq 2 n 1 b 2 /P(n,ni,p ,pi,p2)). 
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Setting m = [8n] , p\ = [/3ip\ , p 2 = [fap] , and p = p - p\ - P2, we obtain 
from (|16p the running time 

'03(l-/3i-A)+(«-l)/8i+(g-2)^a/ 9 
q*I I \/9i+2/3 2 / 



\/3i+/3 2 /\ fe /Vl/ V2^ 

for polynomial size families 3". In particular, we obtain time 0*(3.3432 p ) for 
q = 3 with /3i = 0.281509 and fa = 0.679622, time <3*(7.2562 p ) for q = 4 
with fa = 0.323262 and /3 2 = 0.612790, time O*(15.072 p ) for q = 5 with 
#L = 0.338614 and /3 2 = 0.582673. 

5. A Determinant Sieve for Edge-Coloring 
This section establishes Theorem [6l 

5.1. Tutte's observation revisited. The edges of G can be colored with 
d colors if and only if there exists a set of d — 1 pairwise edge-disjoint perfect 
matchings in G. Indeed, because the graph is d-regular, each color class 
must be a perfect matching. 

Let us now return to Lemma [TTJ We observe that (|14p in effect gives 
us a multivariate generating function for the perfect matchings in G. Our 
strategy is to introduce d — 1 independent copies of this generating function 
and sieve for edge-disjointness. 

5.2. Fingerprinting and identifiability. Let M = (M\, M%, • • • , Mp) be 

an ordered p-tuple of perfect matchings in G. The domain of M is the set 

(25) d(M) = {(e,i) :eeMi}CEx{l,2,..., P }. 

Observe that |d(M)| = pn/2. Let L be a set oipn/2 labels. A labeling of M 
is a mapping A : d{M) — > L. The labeling is bijective if A is a bijection. 

The sieve operates over a multivariate polynomial ring with the coefficient 
field F 2 t and the following indeterminates. Associate with each pair (e, i) G 
Ex{l, 2, ... ,p} an indeterminate x e ^. Associate with each pair (e,l) € ExL 
an indeterminate y e ^. 

The monomial of a labeled p-tuple (M, A) is 

(26) m(M,A)= H x^i£ A(Cji) . 

(e,i)ed(M) 

Lemma 14 (Identifiability). T/ie monomial m(M , A) uniquely determines 
M. Furthermore, if M consists of pairwise edge-disjoint perfect matchings 
and A is bijective, then m(M, A) uniquely determines A. 

5.3. Sieving for bijective labelings. Denote by L the set of all labeled 
p-tuples of perfect matchings of G. For J C L, denote by £[«/] the subset 
of labeled p-tuples of perfect matchings of G whose labeling avoids each 
label in J. Denote by S the set of all bijectively labeled p-tuples of perfect 
matchings of G. 

By the principle of inclusion-exclusion, 

(27) Y, rn(M,\) = J2(-V lJl £ m(M,A). 

(M,A)eS J ^ L (A?,A)e£[J] 
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5.4. Fingerprints of bijectively labeled non-disjoint p-tuples cancel. 

Let 23 be the set of all bijectively labeled p-tuples of perfect matchings of G. 
Partition 23 into 23 = "PU31, where CP is the set of bijectively labeled p-tuples 
of perfect matchings that are pairwise edge-disjoint, and 31 is the set of bi- 
jectively labeled p-tuples of perfect matchings for which there exists at least 
one edge that occurs in at least two matchings in the tuple. Accordingly, 
we have 

Y m(A, h A) = Y m( ^' A ) + Yl m( ^' A ) • 

(M,X)&B (M,A)eIP (M,A)sK 

By a pairing argument essentially identical to the one given in §3.4|, the 
rightmost sum vanishes in characteristic 2. 

5.5. The algorithm. First, the procedure evaluates the following generat- 
ing function for a random assignment of values to the indeterminates. From 
([271) and 3531 we have 

(28) ^Tm(M,A) = ]T(-l)l J l Yl rn(M,X)- 
(M,X)eT J( = L {M,X)eL[J] 

The left-hand side of (|28p is a multivariate polynomial of degree at most 
2pn. It follows from Lemma UM that the polynomial is not identically zero 
if and only if G has a set of p pairwise edge-disjoint perfect matchings. 

It remains to evaluate the right-hand side of (|28|) . Let J C L be fixed. 
The procedure relies on Tutte's Lemma (Lemma II ip . For i = 1,2, ...,p 
define the symmetric n x n matrix T^'(J) as follows. Index the rows and 
columns by the vertices V of G. Define the entries of T^(J) for all u,v 6 V 
by 

(29) e« = {° v if» = »or{«,„}^; 

( x {u,v},i 2^££L\J V{u,v},£ It {U,V} £ E. 

From Lemma [TT1 we have 

v 

(30) Y rn(M,X) = Y[detT^(J). 

(Al,A)G£[J] i=1 

Consequently, for any given assignment of values in F 2 i> to the indetermi- 
nates x £) i and y e ^, the procedure evaluates the left-hand side of (f30|) using a 
total of 0{prfi) arithmetic operations in F 2 &. Taking the sum over J C L, we 
obtain an evaluation of the right-hand side of (j28|) . Taking b = [log 2 4pn] , 
we witness a set of p pairwise edge-disjoint perfect matchings in G as a 
nonzero evaluation of the left-hand side of (|28|) with probability £2(1) in 
time 0*(2 pn ' 2 ) and space polynomial in n. Taking p = d — 1, we obtain 
a polynomial-space randomized algorithm for deciding whether a d-regular 
graph admits a coloring of its edges with d colors in 0*(2^ d ~ 1 ' n ' 2 ) time. 

5.6. Graphs that are not regular. Let m = \E\. We can modify the 
previous algorithm to run in time 0*(2 m ) and space polynomial in n on 
graphs that are not regular. In particular, instead of perfect matchings 
consider matchings, set \L\ = m, in ([29]) set the diagonal entries equal to 1, 
and set p = A, where A is the maximum degree of a vertex in G. 
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